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A model independent expression for the electromagnetic corrections to a phenomenolog- 
ical hadronic pion-nucleon (nN) scattering length a^, extracted from pionic hydrogen, 
is obtained. In a non-relativistic approach and using an extended charge distribution, 
these corrections are derived up to terms of order oP log a in the limit of a short-range 
hadronic interaction. We infer i^-^ = 0.0870(5)m^ which gives for the nNN coupling 
through the GMO relation 5^±p^/(47r) = 14.04(17). 

Keywords: ttN scattering lengths; pionic hydrogen; extended charge distribution and 
short-range hadronic interaction. 

1. Introduction 

Precise Icnowledge of the strong interaction amplitudes at zero energy is important. 
It gives strong constraints to the chiral physics approaclP to QCD at low energy. 
Zero-energy amplitudes can also be used as subtraction constants in dispersion re- 
lations, as for instance in the zero energy pion-nucleon forward dispersion relation^ 
the GMO sum rule, which allows a precise determination of the ttNN coupling^. 
Accurate knowledge of zero energy amplitudes requires careful analysis of electro- 
magnetic corrections of precise experimental information from hadronic atoms. For 
pionic hydrogen recent experimental results give for the strong interaction energy 
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shifPSI 

_ f ( - 7.108 ± 0.013 (stat) ± 0.034 (syst)) eV (precisioiP. ± 0.5%) , , 
^^"^ ~ I ( - 7.120 ± 0.008 (stat) ± °;™g (syst)) eV (precisioiffl ± 0.2%) 

and for the total decay widtlP 

Ti, = [0.868 ± 0.040 (stat) ± 0.038 (syst) ] eV (precisioiP ± 6%). (2) 
It is well knowiJS^ 

eis/Es = iet/EB)il + S^) = Ama a{\ + 5^s)- (3) 

Here Eb = —mo? 12 is the Bohr energy, a = 1/137.036 the fine structure constant, 
TO the 7r~p reduced mass, to = m,^-mpl(mT^- + nip). The correction is the 
deviation from the lowest order shift 

= -(4V2TO)4(0)a. (4) 

The function 0b (r) is the non-relativistic Is Bohr wave function of a point charge. 
The scattering length a corresponds to the elastic threshold scattering amplitude 
defined in the absence of the Coulomb field. 

It is important to understand b\s with an accuracy matching the high experi- 
mental precision. The standard approaclP uses a numerical resolution of coupled 
channel equations. The derivation is model dependent and not fully consistent with 
the ttA^ low-energy expansioiPI. It gives Sis{Sigg) = (—2.1 ± 0.5)%. Calculation 
have also been performed in the framework of QCD+QED with effective field the- 
ory (EFT) techniques. The correction from isospin symmetric pure QCD, evaluated 
within Chiral perturbation Theory (ChPT), is ^= (-4.3 ±2.8)% at leading ordeP 
and (5e = (—7.2 ±2.9)% at next to leading ordeW The large uncertainty in is due 
to the lack of knowledge on one of the low-energy constant entering in the EFT. The 
above determinations consider Sis as corrections to an isospin symmetric world. 

In the present work we study the relation between the strong atomic energy 
shift to the scattering length a'* defined as the one which would be obtained if the 
Coulomb field of the extended charge could be withdrawn and considered as coming 
from an external source. We give here the main features and results of our model, 
more details can be found in Ref. '9. If the Is finite size electromagnetic (e.m.) 
binding energy is denoted by Efs and the total Is binding with strong interaction 
and finite size by E then the strong interaction shift is eis = E — Efg. The non- 
relativistic wave number are kb, Kfs and k for Eb, Efs and E respectively. We 
have — 2mEB and — 2mE. We also use the Bohr radius tb = = (ma)^^ . 

2. Model for the ■7T~p Atom 

In our approach the isospin symmetry is not assumed. The ir^p atom being non- 
relativistic we consider our problem as a non-relativistic quantum one. We work in 
the configuration space. The pion-proton charge distributions are folded to give the 
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Coulomb potential Vc{r). If 14(r) = 0, the low-energy expansion for the angular 
momentum ^ = is given in term of the phase shift Si^o and the momentum q by 

(tan (5fco)/9 = + 9' + - (5) 
Here a'' is the hadronic scattering length and b'^ the range parameter. 

2.1. Toy model 

We first consider a soluble toy model which will be extended to the realistic case. 
We start by studying the case of a single channel. We suppose the charge to be 
concentrated on a sphere of radius R outside the range of the strong interaction. 
Inside this shell, the Coulomb potential is constant with Vc{R) = —a/R ~ 1.4 MeV 
for a i? ^ 1 fm. The Is binding energy Eb = —mo? 12 ~ —3.2 keV is negligible 
compared to the Coulomb field and to the strong interaction inside the charge 
distribution region although its exact value governs the scale of the atom. The inside 
wave function, r < i?, is a standing wave outside the strong interaction region: 

u{r) = N [sin(gc7')/9c + taniJ^^o cos(gcf')/gc] (6) 

where the wave number = 2ma/R — k} ~ 2ma/ R. The external, r < R,ls wave 
function is a Whittakcr function, 

u{r) = (47r)-i/2^e-'='' {2Kr[(l + (1 - A)(l - 7 - log2Kr)] + (1 - A)/A} (7) 

where X = kb/k and 7 = 0.577... is the Euler constant. In Eq. (TJ we have neglected 
terms of order (1 — A)^ ~ {ama^Y — 10~®. Matching the logarithmic derivative of 
the wave function 10 and Q) at i? gives, replacing a by in Eq. ©, 

5is = -2R/rB + 2(a''/rs)[2 - 7 - log(2ami?)] + {2ma/ R){b^ / a^) (8) 

In the absence of strong interaction, the extended charge wave function at ?' = is 

(l>in{0) = (j)B{0){l — R/rB + ■■■) to the present order in a. It is a better e.m. starting 

function than (I)b{0) in Eq. Q. This explains the first term in Eq. 0. The second 

term is a renormalization coming from the external wave function modified at R 

by the hadronic scattering by a factor 1 + 2a^nia[2 — 7 — log(2ami?)]. This factor 

is insensitive to R. The third term is a new effect, it arises from using the correct 

Q 

energy at the point of interaction. It follows from gauge invarianct ■ ^ . 

2.2. Generalization 

Any strong interaction, with the same near threshold hadronic amplitude and with 
a hadronic range inside R, will give the same answer. The difference between the 
true charge distribution and that of the toy model is small and gives a perturbative 
potential. We apply this perturbation to our soluble model to obtain^ an e.m. 
correction which is independent of R to the present order in a. In Eq. |(SJ) i? is 
replaced by {r)em, ^/R by {l/r)em and log(mi?) by (log(TOi?))e„i. This result could 
also have been obtained by a direct calculation. The vacuum polarization correction 
i5"P arises from a long-range potential which modifies the wave function at the origin. 
This model independent correction has been calculated in Ref.llOlto be S^^ = 0.48%. 
It agrees with the earlier numerical value obtained in Ref. |7| 
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2.3. Coupled channel 

In addition to the w^p channel, the n^p atom couples to the ir'^n and 7n channels. 
We describe below the two-channel case with a K-matrix formalism. Let us denote 
the charged tt~p channel by c and the neutral one by o. The complex Coulomb 
threshold amplitude can be written as A'^ ^ = ^ + iqo{K^g)^ where go is the 
momentum in the n'^n channel. The bound state condition leads tc0, 

ei, - tT,j2 = -{47r/2m)<j>l{0)Al, [l + 2AI,{1 + 7)^5] (9) 

The low-energy expansion for the hadronic K-matrix gives i^c.c = Qcc + 
bi^^^ql K,^o = flc.o + ^ioill + 90)72 and Ko,o = a^o,o + &o,o9o- Continuity of the 
wave function matrix and its logarithmic derivative at i?, together with the ex- 
tended charge distribution, lead to 



rB TB \ \ 'rBU,J \T I e^a^ (10) 

= (5<''> + (5^ -I- 

and = ri3/(l + P-^) = (47r/m)4(o)qo [a^^Jl + 5r)]' (11) 

with (5r = (5^''V2 ^6'' ^- {ql + q^bto/C^o-'lo) + 1'^- The Panofsky ratitP', P = 
1.546(9). It is the ratio of the 7r° to 7 production cross sections. 

3. Numerical Results 

The e.m. expectations values in Eq. pO|l are calculated using the folded {tt~ ,p) 
charge distribution obtained from the observed form factors as in Ref. 171 The range 
terms values h^-^ = 6^+„ — — 0.031(9)to~^ and b^-^ = b^+p = — 0.058(9)m~'^ are 
from 'kN scattering datJia. From Eqs. lO and (^UJ, including 5""^ one has, 

ei, = ~(47r/2m)(/)|(0)a;,%(l + 5^'''^ + 5^ + + 5^p). (12) 

Using then the experimental value ei^ of Rcf. 3, plus a two-step iteration from 
a starting value of a^'^, lead to the Sis results given in Table 1. The main source of 
uncertainty comes from the empirical range parameters. The resolution of Eq. (|12|l 



Table 1. Coulomb correction in percent as described in the text. is added to Sia. 





Extended charge 


Renormalization 


Gauge term 


Vacu. Polar. 


Total 
















-0.853(8) 


0.701(4) 


-0.95(29) 


0.48 


-0.62(29) 


St 


-0.427(4) 


0.701(4) 


0.50(23) 


0.24 


1.02(23) 



and that of Eq. (|ll|l . with the experimental value of Fi^ from (O, together with the 
corrections as given in Table 1, plus a negative sign for o^q, lead to 

a;_p = a;?^, = 0.0870(5)m-i ; a^p_,„„ = a^^^^ = -0.0125(4)m-i. (13) 
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The result for a'^-^ is 1.5% smaller and outside the quoted uncertainty of the calcu- 
lated value 0.0883(8)r7i^^ of Ref. |S| using Sigg analysiJ^. Our value of a^-p^^o„ 
is smaller by 2.4% from that of Ref. — 0.128(6)to~^. The ajj_ can be ana- 
lyzed together with a^- jj to give an isovector scattering length (a^r-p — a7r-„)/2. 
We follow Ref. "7 with two small additions: we consider the full triple scattering 
correction-'^ and the additional Fermi motion correction arising from the energy 
dependence of the S-wave isoscalar amplitude at threshold, as first discovered in a 
chiral approacl^l. One obtains (a^-p — a^-„)/-\/2 = — 0.125(l)m~^ in perfect agree- 
ment with the value given in Note here that, in the limit of charge symmetry, 
the further e.m. contributions, for instance from processes such as the 7iV(7A) 
channels and their cross terms, cancel. There is no evidence of isospin violation in 
the isovector amplitude at the current level of accuracy. 

Our result for a^.^, together with the reanalysis of the D scattering length, 
as just described above, allows an improved determination of the ttNN cou- 
pling constant via the GMO sum-rule as given in Eq. (4) of Ref. 2 We obtain 
ff^±pn/(4i') = 14.07(17) to be compared to 14.11(19) in Ref.|2| 

4. Comparison to Previous Approaches and Conclusions 

In former analytical 

approacheP^ii^iei 

using wave functions, the effect of the e.m. 
finite size and the issue of the correct energy of the interaction were basically not 
taken into account. The correction proportional to (a^ is obtained to leading or- 
der a log a, but the part of order a is incorrectly given, leading to a small numerical 
difference. The numerical approach of Ref. 13 with coupled channel Klein-Gordon 
equations, includes correctly the e.m. finite size, the vacuum polarization and the 
renormalization term. It does not have the proper range parameters and makes 
model dependent corrections for isospin violation and radiative effects. The numer- 
ical result for a single channel does not have these problems and agrees with our 
result when using the same input. 

The energy shift ei^ is related to the scattering length in pure QCD in Refs. Q] 
andlHl The calculations are done in a perturbative expansion to the same order in a 
as here in an effective quantum field theory framework. The e.m. EFT corrections 
for mass splitting are implicitly contained in our a^. One cannot then compare 
in detail both expansions. Our renormalization term proportional to (a'*)^ is also 
present in the EFT approach. The main difference is the wave function and energy 
shift corrections linear and inverse in the e.m. charge radius, respectively. No such 
terms seem to be present in the EFT in which only even power (r^")em appear. It 
is due to a different treatment of the e.m. charge form factoJ^. We conjecture that 
the EFT approach should produce terms corresponding to ours as form factors are 
generated using such descriptions. 

Let us conclude by the following important observation. The Coulomb potential 
for the extended charge distribution is regular at short distance and it plays the role 
of an externally applied binding potential in addition to the strong interaction. The 
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problem can then be solved exactly for a model where perturbation can be applied. 
We achieve an accurate determination of the correction terms to the Deser-Trueman 
formula Q in a non-relativistic picture with three quite understood physical effects. 
First the wave function at the origin should match with an extended charge distri- 
bution including vacuum polarization. Second, the correct long-range comportment 
of the wave function changes in a characteristic way the wave function at the ori- 
gin. It leads to a correction proportional to am(a'*)^ log a to leading order. Thirdly 
the low-energy expansion of the scattering amplitude gives a characteristic "gauge" 
correctiorH This effect is as important as the others. Our results should not change 
as long as the hadronic range is within that set by the charge distribution. Our 
approach is general and can be easily applied to other atomic systems like 7r~7r+. 
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